Comparative Investigation of the Microstructure of MgCl2 Aqueous Solutions Using Different X-ray Scattering Sources, Raman Spectroscopy, and Atomistic Simulations.

Aqueous solutions of magnesium chloride (MgCl2(aq)) are often used to test advances in the theory of electrolyte solutions because they are considered an ideal strong 2:1 electrolyte. However, there is evidence that some ion association occurs in these solutions, even at low concentrations. Even a small ion-pairing constant can have a significant impact on the chemical speciation of ions, so it is important to determine whether ion pairing actually occurs. In this study, MgCl2(aq) with concentrations ranging from 1 to 35% was studied using three methods: X-ray scattering (XRS) with the Shanghai Synchrotron Radiation Facility (SSRF) and silver-anode laboratory sources, Raman spectroscopy, and molecular dynamics (MD) simulations with the COMPASS-II and Madrid force fields. XRS results were analyzed in the framework of PDF theory to obtain the reduced structure function F(Q) and the reduced pair distribution function G(r). The F(Q) values from synchrotron radiation and laboratory sources both showed that the tetrahedral hydrogen bonds in bulk water were destroyed with the increased MgCl2 concentration. The results of G(r) indicated that the main peaks centered at 2.05 and 2.80 Å can be ascribed to the interactions of Mg-O and O-O, respectively. The peak at 3.10 Å is attributed to the combined effect of O-O and Cl-O. By comparing the structural information on MgCl2 solution obtained from the two light sources, it was found that both SSRF and silver-anode laboratory sources can reflect the above-mentioned structural information on MgCl2 solution. The radial distribution function (RDF) obtained from MD simulations of MgCl2 solutions assigned the peaks at 2.0, 2.8, and 3.2 Å to the Mg-O, O-O, and Cl-O interatomic pairs, respectively. The decrease in the O-O coordination number confirms that the hydrogen-bonding network of water is disrupted by increasing MgCl2 observed by X-ray scattering. The proportion of Mg-Cl contact ion pairs gradually increases with MgCl2 concentration as does the coordination number. Raman spectroscopy results show that the bond type changes from double donor double acceptor (DDAA) to single donor-single acceptor (DA) with increasing concentration, providing explicit details of the hydrogen-bond evolution in the aqueous solution.

The impact of stoichiometry on the initial steps of crystal formation: Stability and lifetime of charged triple-ion complexes.

Minerals form in natural systems from solutions with varying ratios of their lattice ions, yet non-stoichiometric conditions have generally been overlooked in investigations of new formation (nucleation) of ionic crystals. Here, we investigated the influence of cation:anion ratio in the solution on the initial steps of nucleation by studying positively and negatively charged triple ion complexes and subsequent particle size evolution. Our model systems are carbonates and sulfates of calcium and barium, as it was recently shown that solution stoichiometry affects the timing and rate of their nucleation. Molecular dynamics (MD) simulations and dynamic light scattering (DLS) flow experiments show that nucleation correlates with the stability and lifetime of the initial complexes, which were significantly impacted by the cation:anion stoichiometry and ion type. Specifically, B a S O 4 2 2 - ${{\rm B}{\rm a}{\left({{\rm S}{\rm O}}_{4}\right)}_{2}^{2-}}$ was found to have higher association constants and its lifetime was twofold longer than B a 2 S O 4 2 + ${{{\rm B}{\rm a}}_{2}{{\rm S}{\rm O}}_{4}^{2+}}$ . Similar trends were observed for B a C O 3 ${{{\rm B}{\rm a}{\rm C}{\rm O}}_{3}}$ and C a S O 4 ${{{\rm C}{\rm a}{\rm S}{\rm O}}_{4}}$ . Contrastingly, for C a C O 3 ${{{\rm C}{\rm a}{\rm C}{\rm O}}_{3}}$ , C a C O 3 2 2 - ${{\rm C}{\rm a}{\left({{\rm C}{\rm O}}_{3}\right)}_{2}^{2-}}$ was found to have lower association constants and its lifetime was shorter than C a 2 C O 3 2 + ${{{\rm C}{\rm a}}_{2}{{\rm C}{\rm O}}_{3}^{2+}}$ . These trends in stability and lifetime follow the same asymmetrical behaviour as observed experimentally for particle formation using techniques like DLS. This suggests a causal relationship between the stability and lifetime of the initial charged complexes and the nucleation under non-stoichiometric conditions.

Molecular Insights into the Stability of Titanium in Electrolytes Containing Chlorine and Fluorine Ions.

Titanium and its alloys are protected by a compact and stable passive film, which confers resistance to corrosion by the primary halogen chloride (Cl-) while being less effective against fluoride (F-). Although researchers have recognized different macroscopic corrosion effects of these halide ions on titanium, the underlying mechanisms remain largely unexplored. In this work, the bonding of Cl-/F- with stable passive films was studied in neutral and acidic (pH = 2.3) conditions. The synergistic effect between the interfacial hydrogen bond (HB) structure and halogens on titanium corrosion was first revealed using first-principles calculation and Raman spectroscopy. F- forms more stable halogen-Ti bonds than Cl-, resulting in titanium degradation. The proton combined with F- exhibits a specific synergistic effect, causing corrosion of the passive film. The water hydrogen bond transformation index (HBTI) at the titanium/aqueous interface was 1.88 in an acidic solution containing F-, significantly higher than that in neutral/acid solutions containing Cl- (1.80/1.81) and a neutral solution containing F- (1.81). This work clarifies the structure-activity relationship between HBTI and the destruction of titanium passive films. We propose that the microstructure of the interfacial HB is an undeniable factor in the corrosion of titanium.

Silicon Radical-Induced CH4 Dissociation for Uniform Graphene Coating on Silica Surface.

Due to the manufacturability of highly well-defined structures and wide-range versatility in its microstructure, SiO2 is an attractive template for synthesizing graphene frameworks with the desired pore structure. However, its intrinsic inertness constrains the graphene formation via methane chemical vapor deposition. This work overcomes this challenge by successfully achieving uniform graphene coating on a trimethylsilyl-modified SiO2 (denote TMS-MPS). Remarkably, the onset temperature for graphene growth dropped to 720 °C for the TMS-MPS, as compared to the 885 °C of the pristine SiO2 . This is found to be mainly from the Si radicals formed from the decomposition of the surface TMS groups. Both experimental and computational results suggest a strong catalytic effect of the Si radicals on the CH4 dissociation. The surface engineering of SiO2 templates facilitates the synthesis of high-quality graphene sheets. As a result, the graphene-coated SiO2 composite exhibits a high electrical conductivity of 0.25 S cm-1 . Moreover, the removal of the TMP-MPS template has released a graphene framework that replicates the parental TMS-MPS template on both micro- and nano- scales. This study provides tremendous insights into graphene growth chemistries as well as establishes a promising methodology for synthesizing graphene-based materials with pre-designed microstructures and porosity.

CO2-mineralization and carbonation reactor rig: Design and validation for in situ neutron scattering experiments-Engineering and lessons learned.

CO2 mineralization via aqueous Mg/Ca/Na-carbonate (MgCO3/CaCO3/Na2CO3) formation represents a huge opportunity for the utilization of captured CO2. However, large-scale mineralization is hindered by slow kinetics due to the highly hydrated character of the cations in aqueous solutions (Mg2+ in particular). Reaction conditions can be optimized to accelerate carbonation kinetics, for example, by the inclusion of additives that promote competitive dehydration of Mg2+ and subsequent agglomeration, nucleation, and crystallization. For tracking mineralization and these reaction steps, neutron scattering presents unprecedented advantages over traditional techniques for time-resolved in situ measurements. However, a setup providing continuous solution circulation to ensure reactant system homogeneity for industrially relevant CO2-mineralization is currently not available for use on neutron beamlines. We, therefore, undertook the design, construction, testing and implementation of such a self-contained reactor rig for use on selected neutron beamlines at the ISIS Neutron and Muon Source (Harwell, UK). The design ensured robust attachment via suspension from the covering Tomkinson flange to stabilize the reactor assembly and all fittings (~25 kg), as well as facilitating precise alignment of the entire reactor and sample (test) cell with respect to beam dimension and direction. The assembly successfully accomplished the principal tasks of providing a continuous flow of the reaction mixture (~500 mL) for homogeneity, quantitative control of CO2 flux into the mixture, and temperature and pressure regulation throughout the reaction and measurements. The design is discussed, with emphasis placed on the reactor, including its geometry, components, and all technical specifications. Descriptions of the off-beamline bench tests, safety, and functionality, as well as the installation on beamlines and trial experimental procedure, are provided, together with representative raw neutron scattering results.

Advanced Catalyst Design and Reactor Configuration Upgrade in Electrochemical Carbon Dioxide Conversion.

Electrochemical carbon dioxide reduction reaction (CO2 RR) driven by renewable energy shows great promise in mitigating and potentially reversing the devastating effects of anthropogenic climate change and environmental degradation. The simultaneous synthesis of energy-dense chemicals can meet global energy demand while decoupling emissions from economic growth. However, the development of CO2 RR technology faces challenges in catalyst discovery and device optimization that hinder their industrial implementation. In this contribution, a comprehensive overview of the current state of CO2 RR research is provided, starting with the background and motivation for this technology, followed by the fundamentals and evaluated metrics. Then the underlying design principles of electrocatalysts are discussed, emphasizing their structure-performance correlations and advanced electrochemical assembly cells that can increase CO2 RR selectivity and throughput. Finally, the review looks to the future and identifies opportunities for innovation in mechanism discovery, material screening strategies, and device assemblies to move toward a carbon-neutral society.

Adsorption, activation, and conversion of carbon dioxide on small copper-tin nanoclusters.

Carbon dioxide (CO2) conversion to value-added chemicals is an attractive solution to reduce globally accelerating CO2 emissions. Among the non-precious and abundant metals tested so far, copper (Cu) is one of the best electrocatalysts to convert CO2 into more than thirty different hydrocarbons and alcohols. However, the selectivity for desired products is often too low. We present a computational investigation of the effects of nanostructuring, doping, and support on the activity and selectivity of Cu-Sn catalysts. Density functional theory calculations were conducted to explore the possibility of using small Cu-Sn clusters, Cu4-nSnn (n = 0-4), isolated or supported on graphene and γ-Al2O3, to activate CO2 and convert it to carbon monoxide (CO) and formic acid (HCOOH). First, a detailed analysis of the structure, stability, and electronic properties of Cu4-nSnn clusters and their ability to absorb and activate CO2 was considered. Then, the kinetics of the gas phase CO2 direct dissociation on Cu4-nSnn to generate CO was determined. Finally, the mechanism of electrocatalytic CO2 reduction to CO and HCOOH on Cu4-nSnn, Cu4-nSnn/graphene and Cu4-nSnn/γ-Al2O3 was computed. The selectivity towards the competitive electrochemical hydrogen evolution reaction on these catalysts was also considered. The Cu2Sn2 cluster suppresses the hydrogen evolution reaction and is highly selective towards CO, if unsupported, or HCOOH if supported on graphene. This study demonstrates that the Cu2Sn2 cluster is a potential candidate for the electrocatalytic conversion of the CO2 molecule. Moreover, it identifies insightful structure-property relationships in Cu-based nanocatalysts, highlighting the influence of composition and catalyst support on CO2 activation.

Single-Atom Iridium on Hematite Photoanodes for Solar Water Splitting: Catalyst or Spectator?

Single-atom catalysts (SACs) on hematite photoanodes are efficient cocatalysts to boost photoelectrochemical performance. They feature high atom utilization, remarkable activity, and distinct active sites. However, the specific role of SACs on hematite photoanodes is not fully understood yet: Do SACs behave as a catalytic site or as a spectator? By combining spectroscopic experiments and computer simulations, we demonstrate that single-atom iridium (sIr) catalysts on hematite (α-Fe2O3/sIr) photoanodes act as a true catalyst by trapping holes from hematite and providing active sites for the water oxidation reaction. In situ transient absorption spectroscopy showed a reduced number of holes and shortened hole lifetime in the presence of sIr. This was particularly evident on the second timescale, indicative of fast hole transfer and depletion toward water oxidation. Intensity-modulated photocurrent spectroscopy evidenced a faster hole transfer at the α-Fe2O3/sIr/electrolyte interface compared to that at bare α-Fe2O3. Density functional theory calculations revealed the mechanism for water oxidation using sIr as a catalytic center to be the preferred pathway as it displayed a lower onset potential than the Fe sites. X-ray photoelectron spectroscopy demonstrated that sIr introduced a mid-gap of 4d state, key to the fast hole transfer and hole depletion. These combined results provide new insights into the processes controlling solar water oxidation and the role of SACs in enhancing the catalytic performance of semiconductors in photo-assisted reactions.

The carbon chain growth during the onset of CVD graphene formation on γ-Al2O3 is promoted by unsaturated CH2 ends.

Chemical vapor deposition of methane onto a template of alumina (Al2O3) nanoparticles is a prominent synthetic strategy of graphene meso-sponge, a new class of nano porous carbon materials consisting of single-layer graphene walls. However, the elementary steps controlling the early stages of graphene growth on Al2O3 surfaces are still not well understood. In this study, density functional theory calculations provide insights into the initial stages of graphene growth. We have modelled the mechanism of CH4 dissociation on the (111), (110), (100), and (001) γ-Al2O3 surfaces. Subsequently, we have considered the reaction pathway leading to the formation of a C6 ring. The γ-Al2O3(110) and γ-Al2O3(100) are both active for CH4 dissociation, but the (100) surface has higher catalytic activity towards the carbon growth reaction. The overall mechanism involves the formation of the reactive intermediate CH2* that then can couple to form CnH2n* (n = 2-6) intermediates with unsaturated CH2 ends. The formation of these species, which are not bound to the surface-active sites, promotes the sustained carbon growth in a nearly barrierless process. Also, the short distance between terminal carbon atoms leads to strong interactions, which might lead to the high activity between unsaturated CH2* of the hydrocarbon chain. Analysis of the electron localization and geometries of the carbon chains reveals the formation of C-Al-σ bonds with the chain growing towards the vacuum rather than C-Al-π bonds covering the γ-Al2O3(100) surface. This growth behaviour prevents catalyst poisoning during the initial stage of graphene nucleation.

A Database of Solution Additives Promoting Mg2+ Dehydration and the Onset of MgCO3 Nucleation.

Formed via aqueous carbonation of Mg2+ ions, the crystallization of magnesite (MgCO3) is a promising route to carbon capture and reuse, albeit limited by the slow precipitation of MgCO3. Although magnesite is naturally abundant, forming at low temperature conditions, its industrial production is an energy-intensive process due to the temperatures required to prevent the formation of hydrated phases. The principal difficulty in aqueous conditions arises from the very strong Mg2+···H2O interaction, with high barriers to Mg2+ dehydration. Using atomistic simulations, we have investigated the influence of 30 additive anions (X n-, n = 1-3), ranging from simple halides to more complex molecules, on the first two steps of MgCO3 aggregation from solution, as follows: Mg2+ dehydration and subsequent prenucleative Mg2+···CO3 2- pairing. We have computed the thermodynamic stabilities of solvent shared ion pairs (Mg2+···H2O···X n-) and contact ion pairs (Mg2+···X n-) to reveal the propensity of solution additives to inhibit or promote Mg2+···CO3 2- formation. We have determined the stabilization of undercoordinated hydrated Mg2+ states with a vacant coordination site to which CO3 2- can bind, subsequently initiating MgCO3 nucleation or Mg2+ incorporation into the crystal lattice. Extensive molecular dynamics simulations of electrolyte solutions containing Na2CO3 with different sources of Mg2+ (i.e., MgCl2, MgSO4, and Mg(CH3COO)2) further show that the degree of dehydration of Mg2+ and the structure of prenucleation MgCO3 clusters change depending on the counterion identity. Through a fundamental understanding of the role of solution additives in the mechanism of Mg2+ dehydration, our results help to rationalize previously reported experimental observation of the effect of solvation environments on the growth of magnesite. This understanding may contribute to identifying the solution composition and conditions that could promote the low-temperature CO2 conversion into MgCO3 at industrially relevant scales.

Porous nanographene formation on γ-alumina nanoparticles via transition-metal-free methane activation.

γ-Al2O3 nanoparticles promote pyrolytic carbon deposition of CH4 at temperatures higher than 800 °C to give single-walled nanoporous graphene (NPG) materials without the need for transition metals as reaction centers. To accelerate the development of efficient reactions for NPG synthesis, we have investigated early-stage CH4 activation for NPG formation on γ-Al2O3 nanoparticles via reaction kinetics and surface analysis. The formation of NPG was promoted at oxygen vacancies on (100) surfaces of γ-Al2O3 nanoparticles following surface activation by CH4. The kinetic analysis was well corroborated by a computational study using density functional theory. Surface defects generated as a result of surface activation by CH4 make it kinetically feasible to obtain single-layered NPG, demonstrating the importance of precise control of oxygen vacancies for carbon growth.

Bridging atomistic simulations and thermodynamic hydration models of aqueous electrolyte solutions.

Chemical thermodynamic models of solvent and solute activities predict the equilibrium behavior of aqueous solutions. However, these models are semi-empirical. They represent micro-scale ion and solvent behaviors controlling the macroscopic properties using small numbers of parameters whose values are obtained by fitting to activities and other partial derivatives of the Gibbs energy measured for the bulk solutions. We have conducted atomistic simulations of aqueous electrolyte solutions (MgCl2 and CaCl2) to determine the parameters of thermodynamic hydration models. We have implemented a cooperative hydration model to categorize the water molecules in electrolyte solutions into different subpopulations. The value of the electrolyte-specific parameter, k, was determined from the ion-affected subpopulation with the lowest absolute value of the free energy of removing the water molecule. The other equilibrium constant parameter, K1, associated with the first degree of hydration, was computed from the free energy of hydration of hydrated clusters. The hydration number, h, was determined from a reorientation dynamic analysis of the water subpopulations compared to bulk-like behavior. The reparameterized models [R. H. Stokes and R. H. Robinson, J. Solution Chem. 2, 173 (1973) and Balomenos et al., Fluid Phase Equilib. 243, 29 (2006)] using the computed values of the parameters lead to the osmotic coefficients of MgCl2 solutions that are consistent with measurements. Such an approach removes the dependence on the availability of experimental data and could lead to aqueous thermodynamic models capable of estimating the values of solute and solvent activities as well as thermal and volumetric properties for a wide range of compositions and concentrations.

Optimal Icosahedral Copper-Based Bimetallic Clusters for the Selective Electrocatalytic CO2 Conversion to One Carbon Products.

Electrochemical CO2 reduction reactions can lead to high value-added chemical and materials production while helping decrease anthropogenic CO2 emissions. Copper metal clusters can reduce CO2 to more than thirty different hydrocarbons and oxygenates yet they lack the required selectivity. We present a computational characterization of the role of nano-structuring and alloying in Cu-based catalysts on the activity and selectivity of CO2 reduction to generate the following one-carbon products: carbon monoxide (CO), formic acid (HCOOH), formaldehyde (H2C=O), methanol (CH3OH) and methane (CH4). The structures and energetics were determined for the adsorption, activation, and conversion of CO2 on monometallic and bimetallic (decorated and core@shell) 55-atom Cu-based clusters. The dopant metals considered were Ag, Cd, Pd, Pt, and Zn, located at different coordination sites. The relative binding strength of the intermediates were used to identify the optimal catalyst for the selective CO2 conversion to one-carbon products. It was discovered that single atom Cd or Zn doping is optimal for the conversion of CO2 to CO. The core@shell models with Ag, Pd and Pt provided higher selectivity for formic acid and formaldehyde. The Cu-Pt and Cu-Pd showed lowest overpotential for methane formation.

Sulfate and Molybdate Incorporation at the Calcite-Water Interface: Insights from Ab Initio Molecular Dynamics.

Sulfur and molybdenum trace impurities in speleothems (stalagmites and stalactites) can provide long and continuous records of volcanic activity, which are important for past climatic and environmental reconstructions. However, the chemistry governing the incorporation of the trace element-bearing species into the calcium carbonate phases forming speleothems is not well understood. Our previous work has shown that substitution of tetrahedral oxyanions [XO4]2- (X = S and Mo) replacing [CO3]2- in CaCO3 bulk phases (except perhaps for vaterite) is thermodynamically unfavorable with respect to the formation of competing phases, due to the larger size and different shape of the [XO4]2- tetrahedral anions in comparison with the flat [CO3]2- anions, which implied that most of the incorporation would happen at the surface rather than at the bulk of the mineral. Here, we present an ab initio molecular dynamics study, exploring the incorporation of these impurities at the mineral-water interface. We show that the oxyanion substitution at the aqueous calcite (10.4) surface is clearly favored over bulk incorporation, due to the lower structural strain on the calcium carbonate solid. Incorporation at surface step sites is even more favorable for both oxyanions, thanks to the additional interface space afforded by the surface line defect to accommodate the tetrahedral anion. Differences between sulfate and molybdate substitutions can be mostly explained by the size of the anions. The molybdate oxyanion is more difficult to incorporate in the calcite bulk than the smaller sulfate oxyanion. However, when molybdate is substituted at the surface, the elastic cost is avoided because the oxyanion protrudes out of the surface and gains stability via the interaction with water at the interface, which in balance results in more favorable surface substitution for molybdate than for sulfate. The detailed molecular-level insights provided by our calculations will be useful to understand the chemical basis of S- and Mo-based speleothem records.

Hydrogen-Bond Structure and Low-Frequency Dynamics of Electrolyte Solutions: Hydration Numbers from ab Initio Water Reorientation Dynamics and Dielectric Relaxation Spectroscopy.

We present an atomistic simulation scheme for the determination of the hydration number (h) of aqueous electrolyte solutions based on the calculation of the water dipole reorientation dynamics. In this methodology, the time evolution of an aqueous electrolyte solution generated from ab initio molecular dynamics simulations is used to compute the reorientation time of different water subpopulations. The value of h is determined by considering whether the reorientation time of the water subpopulations is retarded with respect to bulk-like behavior. The application of this computational protocol to magnesium chloride (MgCl2 ) solutions at different concentrations (0.6-2.8 mol kg-1 ) gives h values in excellent agreement with experimental hydration numbers obtained using GHz-to-THz dielectric relaxation spectroscopy. This methodology is attractive because it is based on a well-defined criterion for the definition of hydration number and provides a link with the molecular-level processes responsible for affecting bulk solution behavior. Analysis of the ab initio molecular dynamics trajectories using radial distribution functions, hydrogen bonding statistics, vibrational density of states, water-water hydrogen bonding lifetimes, and water dipole reorientation reveals that MgCl2 has a considerable influence on the hydrogen bond network compared with bulk water. These effects have been assigned to the specific strong Mg-water interaction rather than the Cl-water interaction.

Density functional theory based molecular dynamics study of solution composition effects on the solvation shell of metal ions.

We present an ab initio molecular dynamics study of the alkali metal ions Li+, Na+, K+ and Cs+, and of the alkaline earth metal ions Mg2+ and Ca2+ in both pure water and electrolyte solutions containing the counterions Cl- and SO42-. Simulations were conducted using different density functional theory methods (PBE, BLYP and revPBE), with and without the inclusion of dispersion interactions (-D3). Analysis of the ion-water structure and interaction strength, water exchange between the first and second hydration shell, and hydrogen bond network and low-frequency reorientation dynamics around the metal ions have been used to characterise the influence of solution composition on the ionic solvation shell. Counterions affect the properties of the hydration shell not only when they are directly coordinated to the metal ion, but also when they are at the second coordination shell. Chloride ions reduce the sodium hydration shell and expand the calcium hydration shell by stabilizing under-coordinated hydrated Na(H2O)5+ complexes and over-coordinated Ca(H2O)72+. The same behaviour is observed in CaSO4(aq), where Ca2+ and SO42- form almost exclusively solvent-shared ion pairs. Water exchange between the first and second hydration shell around Ca2+ in CaSO4(aq) is drastically decelerated compared with the simulations of the hydrated metal ion (single Ca2+, no counterions). Velocity autocorrelation function analysis, used to probe the strength of the local ion-water interaction, shows a smoother decay of Mg2+ in MgCl2(aq), which is a clear indication of a looser inter-hexahedral vibration in the presence of chloride ions located in the second coordination shell of Mg2+. The hydrogen bond statistics and orientational dynamics in the ionic solvation shell show that the influence on the water-water network cannot only be ascribed to the specific cation-water interaction, but also to the subtle interplay between the level of hydration of the ions, and the interactions between ions, especially those of opposite charge. As many reactive processes involving solvated metal ions occur in environments that are far from pure water but rich in ions, this computational study shows how the solution composition can result in significant differences in behaviour and function of the ionic solvation shell.

Reconsidering Calcium Dehydration as the Rate-Determining Step in Calcium Mineral Growth.

The dehydration of cations is generally accepted as the rate-limiting step in many processes. Molecular dynamics (MD) can be used to investigate the dynamics of water molecules around cations, and two different methods exist to obtain trajectory-based water dehydration frequencies. Here, these two different post-processing methods (direct method versus survival function) have been implemented to obtain calcium dehydration frequencies from a series of trajectories obtained using a range of accepted force fields. None of the method combinations reproduced the commonly accepted experimental water exchange frequency of 10-8.2 s-1. Instead, our results suggest much faster water dynamics, comparable with more accurate ab initio MD simulations and with experimental values obtained using neutron scattering techniques. We obtained the best agreement using the survival function method to characterize the water dynamics, and we show that different method combinations significantly affect the outcome. Our work strongly suggests that the fast water exchange kinetics around the calcium ions is not rate-limiting for reactions involving dissolved/solvated calcium. Our results further suggest that, for alkali and most of the earth alkali metals, mechanistic rate laws for growth, dissolution, and adsorption, which are based on the principle of rate-limiting cation dehydration, need careful reconsideration.

Solid and Aqueous Speciation of Yttrium in Passive Remediation Systems of Acid Mine Drainage.

Yttrium belongs to the rare earth elements (REEs) together with lanthanides and scandium. REEs are commonly used in modern technologies, and their limited supply has made it necessary to look for new alternative resources. Acid mine drainage (AMD) is a potential resource since it is moderately enriched in REEs. In fact, in passive remediation systems, which are implemented to minimize the environmental impacts of AMD, REEs are mainly retained in basaluminite, an aluminum hydroxysulfate precipitate. In this study, the solid and liquid speciation and the local structure of yttrium are studied in high-sulfate aqueous solutions, basaluminite standards, and samples from remediation columns using synchrotron-based techniques and molecular modeling. Pair distribution function (PDF) analyses and ab initio molecular dynamics density functional theory models of the yttrium sulfate solution show that the YSO4+ ion pair forms a monodentate inner-sphere complex. Extended X-ray absorption fine structure (EXAFS) and PDF analyses show that Y is retained by basaluminite, forming a monodentate inner-sphere surface complex on the aluminum hydroxide surface. EXAFS of the column samples shows that more than 72% of their signal is represented by the signal of basaluminite with which YSO4+ forms an inner-sphere complex. The atomic view of the REE configuration in AMD environments could facilitate a deeper research of REE recovery from waste generated in AMD remediation systems.

The full dynamics of energy relaxation in large organic molecules: from photo-excitation to solvent heating.

In some molecular systems, such as nucleobases, polyenes or the active ingredients of sunscreens, substantial amounts of photo-excitation energy are dissipated on a sub-picosecond time scale, raising questions such as: where does this energy go or among which degrees of freedom it is being distributed at such early times? Here we use transient absorption spectroscopy to track excitation energy dispersing from the optically accessible vibronic subsystem into the remaining vibrational subsystem of the solute and solvent. Monitoring the flow of energy during vibrational redistribution enables quantification of local molecular heating. Subsequent heat dissipation away from the solute molecule is characterized by classical thermodynamics and molecular dynamics simulations. Hence, we present a holistic approach that tracks the internal temperature and vibronic distribution from the act of photo-excitation to the restoration of the global equilibrium. Within this framework internal vibrational redistribution and vibrational cooling are emergent phenomena. We demonstrate the validity of the framework by examining a highly controversial example, carotenoids. We show that correctly accounting for the local temperature unambiguously explains their energetically and temporally congested spectral dynamics without the ad hoc postulation of additional 'dark' states. An immediate further application of this approach would be to monitor the excitation and thermal dynamics of pigment-protein systems.

Prediction of self-assembly of adenosine analogues in solution: a computational approach validated by isothermal titration calorimetry.

The recent discovery of the role of adenosine-analogues as neuroprotectants and cognitive enhancers has sparked interest in these molecules as new therapeutic drugs. Understanding the behavior of these molecules in solution and predicting their ability to self-assemble will accelerate new discoveries. We propose a computational approach based on density functional theory, a polarizable continuum solvation description of the aqueous environment, and an efficient search procedure to probe the potential energy surface, to determine the structure and thermodynamic stability of molecular clusters of adenosine analogues in solution, using caffeine as a model. The method was validated as a tool for the prediction of the impact of small structural variations on self-assembly using paraxanthine. The computational results were supported by isothermal titration calorimetry experiments. The thermodynamic parameters enabled the quantification of the actual percentage of dimer present in solution as a function of concentration. The data suggest that both caffeine and paraxanthine are present at concentrations comparable to the ones found in biological samples.